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Abstract 

A quantum realization of the Relativistic Harmonic Oscillator is realized in terms 
of the spatial variable x and (the minimal canonical representation). The 
eigenstates of the Hamiltonian operator are found (at lower order) by using 
a perturbation expansion in the constant c~ 1 . Unlike the Foldy-Wouthuysen 
transformed version of the relativistic hydrogen atom, conventional perturba- 
tion theory cannot be applied and a perturbation of the scalar product itself is 
required. 



1 Introduction 

The Relativistic Harmonic Oscillator is probably the simplest relativistic system contain- 
ing bound states, yet it exhibits the typical problems of Relativistic Quantum Mechanics. 
Many papers have been devoted to the solution of this relativistc system |], @, 0, §J 
although the first question is probably to define what we understand by a Relativistic 
Harmonic Oscillator. 

In a previous paper || , we adopted an algebraic method for both defining and solving 
such an oscillator equation which started with the Lie operator algebra 



i — p, [E,p] = imu 2 hx, [x,p] = ih(l-\ 

m mc 1 



which is an affine version of the Lie algebra of the 1+1 anti-de Sitter group (~ sl(2, R)) and 
reproduces the appropriate limits (in the sense of Inonii and Wigner group contractions 
0) as w — > (going to the Poincare group in 1+1D) and/or c — > oo (going to the 
Harmonic Oscillator group), although in this paper we will be concerned only with the 

1 Instituto Carlos I de Fi'sica Teorica y Computacional, Facultad de Ciencias, 
Universidad de Granada, Campus de Fuentenueva, Granada, Spain E- 18002 

2 Dipartimento di Scienze Fisiche, Universita degli studi di Napoli, 
Mostra d'Oltremare, Pad. 19, 80125 Napoli, Italy 

3 Instituto de Ffsica Corpuscular, Centro Mixto Universidad de Valencia-CSIC, Burjasot, Valencia, 
Spain E-46100 



1 



c — > oo limit, i.e. the one leading to the Harmonic Oscillator. The solution was given 
through a manifestly covariant representation of the corresponding group (i.e. the wave 
functions are solutions of a Klein-Gordon-like equation). The energy eigenfunctions in 
configuration space consist of a general weight function (the vacuum) which converges to 
a Gaussian in the limit c — > oo , a specific power of the vacuum (reflecting the explicit 
dependence on time of the manifestly covariant representation) which reduces to unity in 
this limit, the usual time-dependent phase factor exp(— iE n t) and a polynomial leading 
to the corresponding non-relativistic Hermite polynomial as c — > oo . As a consequence 
of the explicit dependence on time, the scalar product was defined by using the invariant 
measure dxdt. These representations correspond to a realization in the real variables (x, t) 
of the Discrete series D + (k) of the SL(2, R) group || ||| for Bargmann index k = N = 

In this paper we study the problem appearing in the configuration-space representa- 
tion when the time dependence is factorized out in the search for a minimal, canonical 
realization (given in terms of only x and £) in a way similar to that of the relativistic 
hydrogen atom. In this last case, a Foldy-Wouthuysen transformation leads to a Hamil- 
tonian containing higher-order relativistic corrections to the non-relativistic Hamiltonian 
|T0| . Ordinary perturbation theory provides the energy levels as a power series in c _1 of 
the exact values already obtained from the manifestly covariant (Dirac) equation. 

In the present case, however, the situation is quite different because the terms in 
the Hamiltonian that are considered here as higher-order relativistic corrections to the 
non-relativistic Hamiltonian are not Hermitian with respect to the non-relativistic scalar 
product (with measure dx) and a perturbation of the scalar product itself is required. 
These higher-order relativistic corrections should be understood as a power expansion in 

2 

1/k of the representations of the Discrete series for large Bargmann indices k = N = j^-. 

Our group quantization method |TT| essentially consists in exponentiating the abstract 
algebra (either Poisson or operator algebra) of basic quantities defining a given physical 
system, usually a tt(l)-centrally extended algebra Q, and considering in the resulting group 
(G) the two mutually commuting sets of vectors fields, x R (G) an d X L (G) which generate 
the left and right action of the group on itself. One set of vector fields will constitute 
the operators of the theory, while the other can be used to fully reduce the action of 
the former. To be a bit more precise, we consider the subspace of complex functions on 
G, J-{G), that satisfy the [/(l)-equivariance condition, = i\l/, where H is the central 
generator X^ = X? and ( G U(l). The right-invariant vector fields (generating the 
finite left action) X R e X R {G) ac t 011 ^{G) as a reducible representation (corresponding 
to the Bohr-Sommerfeld quantization). The reduction defining the (true) quantization is 
achieved by imposing on ^F(G) the so-called polarization condition in terms of a subalgebra 
V of left-invariant vector fields, X L ^> = 0,VX L 6 V, which preserve the action of X R 
([X R ,X L ] = 0). The polarization V, originally defined as a maximal left subalgebra of G 
containing the kernel of the Lie algebra cocycle £ and excluding the central generator S, 



llj can be further generalized by allowing it to contain operators in the left enveloping 



algebra U\ L {G) JIJ, [Lj 

The paper is organized as follows: In section 2 the exact solutions for the Relativis- 
tic Harmonic Oscillator, firstly given in [|| and obtained later through a second order 
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polarization in |TJ[], are given. A reduction to a minimal representation (in terms of x 
and 4-) is tempted through a naive elimination of the time variable and the substitution 
ih-jL — > E n — mc 2 , but this leads to a theory which is not unitary. The solution is achieved 
by modifying the scalar product and the operators in a consistent though not quite well 
understood way. In section 3 a perturbative approach to the problem of the reduction is 
proposed. From an exact infinite-order polarization (which represents the "square root" 
of the Casimir leading to the Klein-Gordon equation) we obtain a perturbative expansion 
for a Schrodinger-like equation, where the zeroth-order Hamiltonian is the non-relativistic 
Harmonic Oscillator Hamiltonian. From this Schrodinger equation we inmediately realize 
that the minimal representation is not unitary because the perturbed Hamiltonian (as 
opposite to the case of the hydrogen atom) is not Hermitian with respect to the measure 
dx of the non-perturbed theory. Thus a perturbation of the scalar product itself is pro- 
posed, and the solution coincides (at lower order in powers of 4^) with the exact solution 
proposed in section 2. 



2 The Relativistic Harmonic Oscillator 

To quantize the physical system characterized by the algebra ([!]), we must exponentiate 
this algebra and derive left and right vector fields. The left ones areQ (see [[TJ]] for the 
expression on the right-invariant vector fields): 

~ L p d 2 d P d 

A, = — mu x— — h 



m dx dp mca 2 dt 



~ j Pn d mcx 1 „ 
X — _|_ " 

p mc dp Pq + mc H 

X L = p o d ! P g P mc 1 P 
mc dx mc 2 a 2 dt P + mc h ' 

where Po = V m ° 2 + V 2 + rn 2 uj 2 x 2 and a = y/1 + uj 2 c~ 2 x 2 . 

As was mentioned in the introduction, the representation given by the right-invariant 
vector fields acting on complex wave functions on the group G is reducible. Thus, we 
need to impose conditions on the wave functions (polarization conditions) in order to 
reduce the representation space, and this is achieved by a polarization subalgebra of 
left-invariant vector fields. There is a first-order polarization, leading to the Bargmann- 



Fock representation |]14f , but in this paper we are interested in the configuration-space 



representation, and for this purpose we need a higher-order polarization [|T|, |Tj|. 

There is a second-order polarization which leads to the manifestly covariant represen- 
tation: 

< X?° = (X^) 2 - c 2 m 2 + 2 ^Xt + A^H, Xi > , (3) 

where the numerical parameter A is arbitrary but can be chosen to yield the results pre- 
viously obtained 0. Imposing the £7(l)-equivariance condition and solving the equation 

1 They were called in Ref. M L t ,L p ,L x , respectively 
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Xpty = allows us to factor out the common ( and p-dependence. The remaining equa- 



tion is a Klein-Gordon-like equation for ip(x,t): 

o 2 d f 2 
— 2u x— era 



1 d 2 <p 2imc 2 dtp 2 dtp 2 2 <9 2 y? m 2 c 4 



a 2 dt 2 



ha 2 dt 



dx 



dx 2 h 2 a 2 



<p — A 



mc 2 u 



h 



1 4 

m c 



h 2 



-(p = 



(4) 



In this equation the wave funtion <p(x,t) has the rest mass substracted. If we restore it, 
defining <p = e% mc2t (p, we obtain the more standard expression: 



where □ 



■ dt 2 



2u> 2 x d 
c 2 dx 



□ + 



— a 2 is the D'Alembert operator on anti-de Sitter space- 
time (in 1+1D), R = —2^- is the scalar curvature and x = ^ is a parameter providing 
the coupling of the scalar field to the gravitational field (see [|HJ). 

The normalized (positive-energy) solutions to the Klein-Gordon-like equation are: 



m 2 c 2 



h 1 



X R\<P = 0. 



(5) 



(6) 



where b n = c n = c + n = \ + N x + n = Ejf'V/hw, N x = \ Jl + AN(N - A) 



u ( mu 



T(2N x + l)(2N) n r(N x + ±] 



(7) 



2ttUW 2 n / 2 V^.\ T(2N x + n+l) \J y/~NT(N\) 
N = j^j-, and the Relativistic Hermite polynomials || |i"4"l , H^ N ' X \ satisfy: 

(1 + ^)$H^ - ^(N X + n - \)^H^ + ^(2N X + n)^ = , (8) 

where £ = \J^x (see for the relation between the Relativistic Hermite Polynomials 
and the Gegenbauer polynomials). 

The wave functions are orthonormal according to the t-x scalar product 

<¥ n N ^\¥^> = C™C™ J dxdte-^-^a-^ 2 ^^ x 

H (N,\) H (N,\) g (Q) 

the measure of which comes from the invariant volume Pq dpdxdt after a simple but 
non-trivial regularization of the p- integration ||14|| . 

The annihilation and creation operators (see ||14||), when acting on (p(x,t), have the 
form: 



Z 



h 



2mu 



2mu 



d . ux d mux 
dx c 2 a dt ha 



-iuit 



d ux d mux 

— 1" TT" + 



dx c 2 a dt 



ha 



(10) 
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These operators are the adjoint of each other with respect to the scalar product 
Their action on flBD is: 



£t$(JV,A) 



2N 



n-l 



(n + 1) 



2iV A + n 
2N 



(7V,A) 
n+1 
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The representations here obtained belong to the Discrete series D + {k) of sl(2, R) with 
Bargmann index k = N = Only for half-integer values of k > | they exponentiate 
to a univalued representation of the group SL{2, R) (the rest of the values of k provide, 
however, univalued representations of the universal covering group of SL(2,R)). 

A problem arises, however, when one tries to factorize out the time dependence to 
obtain a minimal representation. The functions a~ Cn H^ N ' x ^ are no longer orthogonal 
unless we modify the scalar product in the form / dx — > J dxa~ 2 (for a discussion on the 
non- fully understood modified scalar product see |14| ). 

With this new scalar product, the normalized wave functions are: 

^(x) = Cj^a-^^(x), (12) 

with 



(J'{N,\) _ / u 



uT fmu\\ 1 r(2iV A + l)(2N) n 
teVhnr) 2 n / 2 V^.\ T{2N x + n + l) \ 



N x + n + 



T(N X 



N x + l \ VNT(N X + 1) 



(13) 

Neither the operators Z and Z^ are adjoint to each other with respect to the new scalar 
product, so they must be appropiately corrected. The relation between the corrected 
operators and the old ones, when acting on the energy eigenfunctions, is given by: 



Z' 



-iujt 



N x + n 
\l N x + n + \ 



l 

5 Z 



Aut 



N x + n+J ^ 



\ N x + 



n 



(14) 



These expresions are the generalization to arbitrary A of the ones given in [14| for A = 1. 
These results suggest that there is a unitary transformation U relating the manifestly 



covariant and the minimal representatations, i.e. 
formation is of the form: 



^W)( x ) = U^ N 'V(x,t). This trans- 



U 



MX) 
\l hujN x 



(15) 



when acting on energy eigenfunctions. Having into account that the Hamiltonian H = 

EV*'*) *W A J(ar,t), we can write U as: 



ih§- t satisfies H ^> x) (x, t) 



U 



H 



\ TiuN x 



(16) 
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when acting on an arbitrary function. With this expresion we can obtain the form of the 
operators Z' and Z'^ when acting on arbitrary functions, not only energy eigenf unctions, 
simply transforming them by If: Z' — UZU^ 1 , Z'* = UZ^U^ 1 . 

One of the problems with this approach is the lack of a Schrodinder-like equation 
providing an expresion of H in terms of x and since we have only at our disposal the 
Klein-Gordon-like equation (|4]) and we would need its "square root" . In the next section 
we shall obtain an expression for H (at low order in 1/N) in terms of x and -r- through 
a Schrodinger-like equation derived from an infinite-order polarization. 

3 A Perturbative Approach Involving a Perturbed 
Scalar Product 

Another way of approaching the t-factorization problem consists in taking the "square 
root" of the second-order polarization above, a solution to the conditions defining an 
(infinite-order) polarization given in a power series < X^°,X^ >, where 



X t °° = X t L + -{ x /ti'u 2 N{N - A) - Ti 2 c 2 



mc 



(17) 



This infinite-order polarization can be solved order by order to obtain a perturbative 
expansion for the wave functions. The first-order terms in 1/N, or, equivalently, in 1/c 2 , 
for X£° are (we are taking into account the other polarization equation, X^^f = 0): 



X? 



X 



0/_(vL\2 ^_ 

2m [ x) AN 



{X L x f-u{l-2a) 



2m 2 u 



Here we have introduced A = % since, although the parameter A can take any value, to 
obtain the correct energy for the non-relativistic Harmonic Oscillator in the limit c — » oo, 
it has to be of order lower or equal to 1/N. In particular, the solutions obtained in ||, 
characterized by A = 1, do not satisfy this requirement since the energy eigenvalues in 
the limit c — > oo are E n = hum, losing the vacuum energy ^huj which characterizes the 
quantum fluctuations of the Harmonic Oscillator system. 

Once the common (- and p-dependences have been factorized out, the new polarization 
gives for ip(x, t): 



h 4 <9V 



8m 3 c 2 dx 4 2m 



1 99/-1 UJ X . 

-mu x (l-—)<p + 



3tu x d if h uj l x dip 
~2& 

h 2 u 2 



+ 



Amc 2 



dx 2 2mc 2 dx 
1 - 2a)ip + 0(c- 4 ) = (H-mc 2 )(p. (19) 



(Note that we have substracted the rest mass from the Hamiltonian, in order to get 
the correct non-relativistic limit). Ordinary perturbation theory, when applied to the 
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Hamiltonian ([19]), yields the correct energy E n « %uj{\ + n) + M(i _ 4a) + 0(c- 4 ) and 
eigenfunctions (p n (x,t) = exp(—jrE n t)<j) n (x), where 

7 , . { muj\ i exp(— mujx 2 l2Ti) f ... 1 r 1 T ... ... 

Wx) = Uf) W J ( g -tt) + g]ylg g ^tt) + g ^K)+ 

4n(n - l)/r„_ 2 (|) - 2n(n - l)(n - 2)(n - 3)/f„_ 4 (f)] + 0(<T 4 )) , (20) 

for any value of a. The differences among different values of a come at order 1/c 4 , as far 
as a is of order 1. If a was of order c 2 (i.e., A of order 1), then the corrective terms would 
appear at order 1 in the energy and at order 1/c 2 in the wave function, as is the case 
of the solutions considered in ||, as commented before. In this paper we shall restrict 
ourselves to the case a of order 1 (A of order 1/c 2 ), since it reproduces completely the 
non-relativistic limit, not only the wave functions, but also the energy. 

However, / dx(p n (x)(j) rn (x) ^ for m ^ n, the reason being that the perturbative terms 
added to the non-relativistic Hamiltonian ([19]) are not Hermitian. It is important to stress 
that with this perturbative approach one can see clearly why the energy wave funtions 
are not orthogonal to each other, since we have an explicit expression for the Hamiltonian 
(although at lower order in powers of 1/c 2 ), and this proves to be non-Hermitian. 

We could add corrective terms to the perturbed Hamiltonian in order to make it 
Hermitian (note that the Hermiticity of H does not depend on the value of a), but 
since the expression of H is given a priori, from a more general theory 2 , we are forced 
to associate the problem with the scalar product, which is no longer adequate for the 
perturbed theory. Thus, a perturbation of the scalar product itself is required, and this 
can be achieved by considering an arbitrary power expansion in c~ 2 as the perturbed 
measure, and determining the corresponding coefficients by using the condition that the 
Hamiltonian be Hermitian at each order. As a result, we get the expression 

dx{l -wVc- 2 + 0(c- 4 )} (21) 

for the perturbed measure, restoring in this way the unitarity of the theory. The power 
series in c~ 2 we have obtained corresponds to that of / dxa~ 2 , in agreement with the 
solution proposed for the exact case in the previous section (see also Jl4| for a discussion). 
The normalized wave functions (according to the perturbed measure) are: 

*„(•'■) = ( 1+ 2n ^ 1 ) '■>»<•'•> 

1 

(x) + 



16N 



^{n + 1) (n + 2) (n + 3) (n + 4)0 n+4 (x) 



+A^{n + l){n + 2)<p n+2 {x) + 8(n + ^n(x) + ^n{n - l)0„_ 2 (x) (22) 



2 The perturbed Hamiltonian has been obtained from a power series expansion for a higher-order 
polarization, but it could well had been obtained from a Foldy-Wouthuysen transformation of a Dirac 
equation for the Relativistic Harmonic Oscillator, as it is the case for the corrective relativistic terms of 
the Hydrogen atom system . 
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-Jn{n - l)(n - 2){n - 3)0 n _ 4 (x) 



where 4> n (x) = (jjrjf) 4 exp ^^^/ 2 ^ i? n (£ ) are the normalized (non-relativistic) Harmonic 
Oscillator wave functions. It is easy to check that these wave functions coincide, up 
to 0(l/c 4 ), with the exact wave functions ( |T2"D normalized with respect to the measure 
/ dx/a 2 . 

We can easily check that the relation between Q(x) and (p(x,t), 

M^ = e^ t (l+ T ^\ Vn (x,t), (23) 

corresponds with the lower order terms of the transformation U given in for energy 
eigenfunctions. But now, when passing to arbitrary functions in (|16|), we have an explicit 
expression, at least at lower order, for the Hamiltonian H in terms of x and j- given by 

& 

We think that the present perturbative approach provides strong support to the idea 
that Perturbation Theory must be in general modified in order to incorporate "non- 
Hermitian" corrections, like terms of the form x 2 -^ } which appear as a consecuence of the 
non-trivial curvature of the space-time (such as the Anti-de Sitter universe). Furthermore, 
our particular example also constitutes an important step towards the solution of the 



general Cauchy problem in non-globally hyperbolic spacetimes [16, T7J, the Anti-de Sitter 
universe being one of the few exactly solvable cases. 
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